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Abstract 

We prove that a model atom having one bound state will be fully 
ionized by a time periodic potential of arbitrary strength r and fre- 
quency oj. Starting with the system in the bound state, the survival 
probability is for small r given by e^^* for times of order ~ r~^", 
where n is the minimum number of "photons" required for ionization 
(with large modifications at resonances). For late times the decay is 
like with the power law modulated by oscillations. As r increases 
the time over which there is exponential decay becomes shorter and 
the power law behavior starts earlier. Results are for a parametrically 
excited one dimensional system with zero range potential but compar- 
ison with analytical works and with experiments indicate that many 
features are general. 

PACS: 32.80 Rm, 03.65 Bz, 32.80 Wr. 
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1 Introduction 



The solution of the Schrodinger equation with a time dependent potential 
leading to transitions between bound and free states of a quantum system 
is clearly of great theoretical and practical interest. Fermi's golden rule 
(based on a deep physical interpretation of first order perturbation the- 
ory) gives the decay exponent of the survival probability for a system in a 
bound state subjected to a weak external oscillating potential, with frequency 
uj > uq = —Ub/h, Ub the energy of the bound state [T]. The approaches used 
to go beyond the golden rule include higher order perturbation theory, semi- 
classical phase-space analysis, Floquet theory, complex dilation, exact results 
for small fields, and numerical integration of the time dependent Schrodinger 
equation 0-^1]. These works have yielded both theoretical understanding 
and good agreement with dissociation experiments in strong laser fields. In 
particular they have been very successful in elucidating much of the rich 
structure found in the experiments on the multiphoton ionization of Ryd- 
berg atoms by microwave fields 0^10] • Explicit results for realistic systems 
require, of course, the use of some approximations whose reliability is not 
easy to establish a priori. 

It would clearly be desirable to have examples for which one could com- 
pute the time evolution of an initially bound state and thus the ionization 
probability for all values of the frequency and strength of the oscillating 
potential to as high accuracy as desired without any uncontrolled approxi- 
mations. This is the motivation for the present work which describes new 
exact results relating to ionization of a very simple model atom by an oscil- 
lating field (potential) of arbitrary strength and frequency. While our results 
hold for arbitrary strength perturbations, the predictions are particularly 
explicit and sharp in the case where the strength of the oscillating field is 
small relative to the binding potential — a situation commonly encountered 
in practice. Going beyond perturbation theory we rigorously prove the exis- 
tence of a well defined exponential decay regime which is followed, for late 
times when the survival probability is already very low, by a power law de- 
cay. This is true no matter how small the frequency. The times required 
for ionization are however very dependent on the perturbing frequency. For 
a harmonic perturbation with frequency u the logarithm of the ionization 
time grows like r~^", where r is the normalized strength of the perturbation 
and n is the number of "photons" required for ionization. This is consistent 
with conclusions drawn from perturbation theory and other methods (the 
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approach in jHj being the closest to ours), but is, as far as we know, the first 
exact result in this direction. We also obtain, via controlled schemes, such 
as continued fractions and convergent series expansions, results for strong 
perturbing potentials. 

Quite surprisingly our results reproduce many features of the experimen- 
tal curves for the multiphoton ionization of excited hydrogen atoms by a 
microwave field These features include both the general dependence of 
the ionization probabilities on field strength as well as the increase in the life 
time of the bound state when —nhu, n integer, is very close to the binding 
energy. Such "resonance stabilization" is a striking feature of the Rydberg 
level ionization curves [3]. These successes and comparisons with analytical 
results suggest that the simple model we shall now describe contains 

many of the essential ingredients of the ionization process in real systems. 

1.1 Description of the model 

The unperturbed Hamiltonian in our model is, 

Hq has a single bound state Ub{y) = i/Po6~p"I^I, po = with energy —htuo = 
—h'^pl/2m and a continuous uniform spectrum on the positive real line, with 
generalized eigenf unctions 

u(k, y) = ( e'^y ^^e'l'^^l^ , -oo < A; < oo 

V27r V Po + Ak\ J 



and energies U^k'^ /2m. 

Beginning at some initial time, say t = 0, we apply a parametric perturb- 
ing potential —gri(t)6{y), i.e. we change the parameter g in Hq to ^'(1 + ri(t)) 
and solve the time dependent Schrodinger equation for ip{y,t), 



^{y, t) = eit)ub{y)e''^'' + / e{k, t)u{k, y)e"'^'dk {t > 0) (2) 

J —oo 

with initial values 6{0) = 1, Q{k, 0) = 0. This gives the survival probability 
as well as the fraction of ejected electrons \Q{k,t)\'^dk with (quasi-) 
momentum in the interval dk. 
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This problem can be reduced to the solution of an integral equation ^7] . 
Setting 



e{t) = l + 2i [ Y{s)ds (3) 
Jo 



t 

2^ 



e(A;,t) = ^ / F(s)e^(i+^')^rfs (4) 



^(1-^1^1) Jo 
Y{t) satisfies the integral equation 

Y{t) = r]{t) |l + j\2t + M{t - t')]Y{t')dt''^ (5) 

where we have set po = ujq = h = 2m = 2 ~ ^ ^^'^ 

, , 2i ^2g-*s(i+«2) I fi e"^" , 

M(s) = — ^ — du=-\- -^du. 

^ ' n Jo l + u^ 2 V Try, u^^ 

2 Results 

Our first exact result is the following: When ri{t) is a trigonometric polyno- 
mial, 

n 

Tfit) = ^[Aj sm{jujt) + Bj cos(ju;t)], (6) 
i=i 

the survival probability \0{t)\'^ tends to zero as t 00, for all u > 0. That 
is there will be full ionization for arbitrary strength and frequency of the 
oscillating field. 

Since the main features of the argument are already present in the sim- 
plest case rj = r sin(ci;t) we now specialize to this case. The asymptotics of 
Y{t) are obtained from its Laplace transform y{p) = e^^^Y{t)dt, which 
satisfies the functional equation (cf. (0)) 

^^^^ 2 \ ^l-tp + cO - 1 y/l-tp-UJ - 1 J + p2 ^ > 

with the boundary condition y{p) — > as Q{p) —>■ ±00 (the relevant branch 
of the square root is {1 — ip — uy/"^ = —i{uj — 1 + ipY^"^ for u > 1). We 
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show that the solution of ((Tj) with the given boundary conditions is unique 
and analytic for 3?(p) > 0, and its only singularities on the imaginary axis 
are square-root branch points. This in turn implies that does indeed 

decay in an integrable way. 

The ideas of the proof carry through directly to the more general periodic 
potential © and we have obtained analogous results for a two delta functions 
reference potential 

Full ionization is in fact expected (for entropic reasons) to hold generically, 
but has, as far as we know, only been proven before for small amplitude of 
the oscillating potential with a; > 1 (IID], [H])) or for sufficiently random 
perturbations f[TU]). 

The detailed behavior of the system as a function of t, u, and r is obtained 
from a precise study of the singularities of y{p) in the whole complex p-plane. 
Here we discuss them for small r; below, the symbol e describes error bounds 
of order 0(r^~^), | < 5 < 1. At p = {inu — i : n E Z}, y has square root 
branch points and y is analytic in the right half plane and also in an open 
neighborhood A/" of the imaginary axis with cuts through the branch points. 
As \^{p)\ ^ oo in ATwe have \y{p)\ = 0{ru;\p\~'^). If - ^\ > C„e, n G Z+, 
for some constants C„, then for small r the function y has a unique pole Pm 
in each of the strips —mu > '^{p) + 1 ± e > —mu — u, m G Z. ^{pm) is 
strictly independent of m and gives the exponential decay of 6. 

The analytic structure of y is indicated in Figure 1 where the dotted lines 
represent (the square root) branch cuts and the dark circles are simple poles. 
The function Y is the inverse Laplace transform of y 



where the contour of integration C can be initially taken to be the imaginary 
axis iR, since y is continuous there and decays like for large p. 

We then show that C can be pushed through the poles, collecting the 
appropriate residues, and along the branch cuts as shown in Figure 1. The 
residue at the pole p^ is proportional to e^^""'"*™'^)* while the (rapidly conver- 
gent) integral along the m— th branch cut is (as seen by standard Laplace 
integral techniques), a function whose large t behavior is Kme^^'^^t~^^'^ (the 
—3/2 power comes from the fact that y has square root branch points; Km 
is some constant). 

A detailed analysis along these lines yields ([IE], jT^)- 
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Figure 1: Singularities of y and relevant inverse Laplace contours. 



oo 

e{t) = e-^("'")*F^(t) + e(™"-')*/i^(t) (9) 

m=— oo 

where is periodic of period 27ra;~^ and hm{t) ~ Ylf=o Cm,j't~^^'^~-' for large t 
in more than a half plane centered on the positive real half-line. Not too close 
to resonances, i.e. when \uj — n^^\ > e, for all integer n, \F^{t) \ = 1 + O(r^) 
and its Fourier coefficients decay faster than rl^'"l|m|~'™''/^. Also, the sum in 
(Q does not exceed 0(r^t~^/^) for large t, and the decrease with m faster 
than r'™'!. 

By Q, for times of order l/F where F = 23fJ(7), the survival probability 
for u not close to a resonance decays as exp(— Ft). This is illustrated in 
Figure 2 where it is seen that for the exponential decay holds up to 

times at which the survival probability is extremely small. Note also the slow 
decay for u = .8, when ionization requires the absorption of two photons. 
For even larger r one can note in the figure oscillatory behavior. This is 
expected from equation ©. 

When r is larger (inset in Fig. 2) the ripples of |-Fa;(^)| are visible and 
the polynomial-oscillatory behavior starts sooner. Since the amplitude of 
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Figure 2: Plot of log^^Q vs. time in units of lUq'^ for several values of to 

and r. The main graph was calculated from (8) and the inset used numerical 
integration of (5). 
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the late asymptotic terms is O(r^) for small r, increased r can yield a higher 
late time survival probability. This phenomenon, sometimes referred to as 
"adiabatic stabilization" [201, associated with the perturbation- 

induced probability of back-transitions to the well. 

Using continued fractions F can be calculated convergently for any uj and 
r. For small r we have 



Ijj 



■IJJ 



F ~ < 



)-2n+2 



\fnjj 



„2n 



- Vl - ma;) 



2 nu; ■ 



if Co" G (1 + £, cxd) 



if G (i + £, 1 -£) 



if u; e (- -he, 



(10) 



m<n 



Eq. (fTIH) agrees with results of perturbation theory. Thus, for > 1, F is 
given by Fermi's golden rule since the transition matrix element between 
the bound state with energy —1 and the continuum state with energy /c^ is 



< Ui,{y)\b{y)u{k,y) > 



1 



27r 1 + F 

while the density of states is 

The behavior of F is different at the resonances uj~^ G N. For instance, 
whereas if u is not close to 1, the scaling of F implied by (fTUI) is when 

when UJ — 1 = r^/ \/2 we find 



a; > 1 and when ^ < a; < 1 



2i/4 2^/4 



16 



In Figure 3 we plot the behavior of F"^, as a function of u, for a small 
value of r. The curve is made up of smooth (roughly self-similar) pieces for u 
in the intervals (n~^, (n — 1)^^) corresponding to ionization by n photons. At 
resonances (for small r these occur for uj~^ close to an integer), the coefficient 
of r^", the leading term in F, goes to zero. This yields an enhanced stability 
of the bound state against ionization by perturbations with such frequencies. 
The origin of this behavior is, in our d = 1 model, the vanishing of the matrix 
element in (fTTj) at k = 0. This behavior should hold quite generally in c? = 3 
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Figure 3: log^^QF^^ as a function of uj/ujq at r = 0.01. At r = 0.2 (inset) 
shift of the resonance is visible. In the inset the continued fraction was used, 
while for the main figure the data was extrapolated from (jlO|) . 



where there is a factor k inT coming from the energy density of states near 
k = 0. As r increases these resonances shift in the direction of increased 
frequency. For small r and uj = 1 the shift in the position of the resonance, 
sometimes called the dynamic Stark effect [1], is about 

In Figure 4 we plot the strength of the perturbation r, required to make 
|6'(t)p = I for a time of about 700 oscillations of the perturbing field, i.e. time 
is measured in units of uj~^, as a function of u. Also included in this figure are 
experimental results taken from Table 1 in [HI, see also Figures 13 and 18 there 
for the ionization of a hydrogen atom by a microwave field for approximately 
the same number of oscillations. In these still ongoing beautiful series of 
experiments, the atom is initially in an excited state with principal 

quantum number Uq ranging from 32 to 90. The "natural frequency" c^o 
is there taken to be that of an average transition from uq to uq ± 1, so 
ti^o ~ '^o^- The strength of the microwave field F is then normalized to the 
strength of the nuclear field in the initial state, which scales like tzq ^. The plot 
there is thus of UqF vs. HqU. To compare the results of our model with the 
experimental ones we had to relate r to n^F, and given the difference between 
the hydrogen atom perturbed by a polarized electric field Vi = xF sm{ujt), 
and our model, this is clearly not something that can be done in a unique 
way. We therefore simply tried to find a correspondence between UqF and r 



9 



0.36 



- 0.12 



0.30 



0.1 



0.24 



0.16 



-Experiment [3] 



0.08 

:v 

0.06 



0.12 



Eq. (8)-i 



0.04 



0.06 



0.02 



0.2 0.4 0.6 



1 CD/ CI),, 



Figure 4: Threshold amphtudes for 50% ionization vs. cj/cuq, calculated from 
eq. (jlOj) . and in the experiment [3]. In the calculation, the dynamic Stark 
effect was approximated using an averaged r over the range. 



which would give the best visual fit. Somewhat to our surprise these fits for 
different values oi uj/ujq all turned out to have values of r close to Srig-F. 

A correspondence of the same order of magnitude is obtained by com- 
paring the perturbation-induced shifts of bound state energies in our model 
and in hydrogen. We note that the maximal values of r in Figure 4 are still 
within the regime where only a few terms in © are sufficient. 

In Figure 5 we plot vs. r for a fixed t and two different values of 

ui. These frequencies are chosen to correspond to the values of oj/uq in the 
experimental curves, Figure 1 in [H] and Figure lb in The agreement is 
very good for uo/ujq .1116 and reasonable for the larger ratio. Our model 
essentially predicts that when the fields are not too strong, the experimental 
survival curves for a fixed n^uj (away from the resonances) should behave 

essentially like exp ^— C[nQF]"o^ ujt^ with C depending on n^oj but, to first 

approximation, independent of n^F. 
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Figure 5: Ionization fraction at fixed t (corresponding to 300 oscillations) as 
a function of amplitude. 

3 Concluding remarks 

Given the simplicity of our model, the similarity (using a minimal numbers 
of adjustable parameters) with the experiments on hydrogen atoms was sur- 
prising to us. As already noted these experimental results and in particular 
the resonances can be understood quite well, including some details, by do- 
ing calculations on the full hydrogen atom or a one dimensional version of it 
^-|6j. Still, it is interesting to see that similar structures arise also in very 
simple models. It suggests that various features of the ionization process 
have a certain universal character. To really pin down the reason for this 
universality will require much further work. 

We note that for u) > ujq, in the limit of small amplitudes r, a predomi- 
nantly exponential decay of the survival probability followed by a power-law 
decay was proved in ^T] for three dimensional models with quite general 
local binding potentials having one bound state, perturbed by r cos{u!t)V (y) , 
where is a local potential. Our results for general u and r can be thought 
of as coming from a rigorous Borel summation of the formal (exponential) 
asymptotic expansion of Y for t — 00. These methods can be extended to 
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other systems [T7] including, we hope, reahstic ones. 
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